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1 Introduction The Aharonov-Bohm (AB) effect shows the importance of the electromagnetic gauge 
field in quantum mechanics and the fact that a single electron has a fundamental unit of magnetic flux 
$0 = 27r/e (units for which h ~ c = 1 are used here). One direct consequence of the AB effect in 
solid state physics is persistent cuiTents [1, 2, 3]. Persistent cuiTents are driven by an external gauge field 
produced when a magnetic flux penetrates a hollow area of a system. The external gauge field breaks the 
time-reversal symmetry and can induce equilibrium persistent currents in the ground state. In a previous 
study [4], a twist-induced gauge field was proposed that also affects persistent currents but still preserves 
the time-reversal symmetry. This gauge field is useful in analyzing persistent cuiTents in toroidal carbon 
nanotubes as well as in other systems [5]. On the other hand, the effect of lattice deformation in carbon 
nanotubes on conducting electrons was studied, and it was found that electrons near the Fermi level are 
affected by a gauge field produced by the deformation [7] (called the "deformation-induced gauge field"). 
Because a deformation depends on position, the deformation-induced gauge field is also position depen- 
dent. The purpose of this paper is to show that the twist-induced gauge field corresponds to a constant part 
of the deformation-induced gauge field, meaning that a local extension of the twist-induced gauge field 
appears naturally in the low energy dynamics of conducting electrons in a deformed carbon nanotube. The 
gauge field presented here provides a unified way of examining the effect of local lattice deformation in 
carbon nanotubes on the local energy gap. At the conclusion of this paper, some problems associated with 
the above-mentioned gauge field are discussed. 

2 Basics of Carbon Nanotubes In this section, the basic ideas used to describe the quantum mechan- 
ical behavior of conducting electrons in carbon nanotubes [6, 7] are explained. These ideas underlie the 
work presented in subsequent sections. 

A carbon nanotube is a graphene sheet folded into a cylinder, and graphene consists of many hexagons 
(in a honeycomb lattice) whose vertices are occupied by carbon atoms. Each carbon atom supplies one 
conducting electron that determines the electrical property of the carbon nanotube. The graphene lattice 
structure is illustrated in Fig. 1 . There are two symmetry translational vectors (ai and 02 satisfying | ai p = 
|a2p — 2ai ■ 02) on the planar honeycomb lattice. The lattice structure of a carbon nanotube can be 
specified by two mutually perpendicular vectors, the chiral Ch and the translational T, defined respectively 
by Ch = nai + ma2 and T = pai + qa2- Here, Ch fixes the lattice structure around the axis and is denoted 
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simply by Ch = (n, m). The corresponding wave vector k can be decomposed as A; = niki + /U2fc2, 
where and ^2 are integer coefficients of vectors ki and k2 satisfying the periodic boundary conditions 
Ch • k\ = 27r, Ch ■ k2 = ^, T ■ ki — Q and T ■ k2 = 27r. Thus, Ch and T specify the kinematics of the 
electron. 

A nontrivial property of a single electron wave function in a honeycomb lattice is that the wave function 
has two components. This is because two sites exist in the unit cell, depicted as circles that are black 
(A-site) or white (B-site) in Fig. 1. The two-component wave function is written as 



I*) 



(1) 



It is noted that this two-component nature does not relate to electron spin. However, there is a (theoretical) 
similarity between them such that the two-component nature is sometimes referred to as a "pseudo-spin." 
For example, when an electron takes a turn (around the Fermi point) in k-space through the impurity 
scattering, the wave function of the electron takes the (Berry) phase shift of tt resulting in an absence of 
back scattering [8]. The sign change of the wave function is similar to that of spin rotation in space. 



Unit Cell 




Fig. 1 Structure of a honeycomb lattice having the 
two symmetry translation vectors ai = v^SaccCa! 
and 02 = {\/?>/2)accex + (3/2)accej,, where 
and Cy are unit vectors in the direction of x and y, 
and Occ is the nearest-neighbor bond length. The 
black (white) circles indicate the A (B) sublattices. 
Vectors Ra (a = 1,2,3) point to the nearest- 
neighbor sites of an A-site. They are given by 



Ri — flccEy, R2 — —{\/Z/2)accex 
and R:, = (^/2)a^^ex - (l/2)ac, 



(l/2)ao 



The Hamiltonian of this system {dynamics) is thus given as a 2 x 2 Hermite matrix. We consider a 
general form of the Hamiltonian: 



(2) 



where TLaa is the probability amplitude from an A-site to the same site, and Hab is that from an A-site to 
the nearest B-site. We assume the nearest-neighbor tight-binding Hamiltonian Ho, with hopping integral 
given by 



0=1,2,3 ieA 



(3) 



where A denotes an A-sublattice, and a\ are canonical annihilation-creation operators of the electron 
at site i satisfying the anticommutation relation {ai, a]} = <^ij, and site i + a indicates nearest-neighbor 
sites (a = 1, 2, 3) of site i. Here, T-Laa and TLbb are neglected because they are assumed to be equivalent 
and to have a constant value, meaning that they do not affect the energy spectrum but only shift the origin 
of the energy. By means of the Bloch theorem, we obtain T-Lab = Hqa = ^ J2a e**^'^". Here, Ra 
(a = 1, 2, 3) are vectors pointing to the nearest neighbor of an A-site (see Fig. 1). By diagonaUzing the 



1,2,3 ' 



, where 



Hamiltonian, we obtain the energy eigenvalue as a function of k as E{k) = ±14- 
the negative sign indicates the valence energy band, and the positive sign denotes the conduction band. It 
is crucial to know that there are two positions in k-space where E{k) vanishes. They are located at ±K 
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where K = {2hi + 62)/3 and called the K and K' points. Here, hi and 62 are reciprocal lattice vectors of 
the graphene defined by • bj — 27r5jj (i, j = 1,2). 

In the following, metallic or semiconducting behavior dependent on the lattice structure around axis Ch 
as exhibited by the carbon nanotube is shown. Assuming a nanotube with long length, that is |r| ^ \Ch\, 
then lines can be drawn in k-space indicating that the energy eigenstates are located on the lines (see 
Fig. 2). The energy gap of the system is determined by the configuration of the line located nearest to 
the Fermi points. Here, the electrons in the lines (the energy bands of the lines) are called "low-energy 
electrons" (low energy modes). When the line cross the K-point, the system exhibits metallic properties 
and otherwise semiconducting properties. This is because the conduction band and valence band touch 
only at the K and K' points. It is important to mention that low energy electrons around the Fermi points 
generally rotate around the axis as opposed to a naive consideration. It is incorrect to think that low energy 
electrons should always have non-rotating modes or that the wave function around the axis is constant. 
The motion of low energy electrons around the axis depends on the lattice structure around the axis. One 
can verify the relationship between the chiral vector and the motion of low energy electrons in Table I of 
reference [4]. 



(a) 



K2 K2 



K' ' 



fcl 



A d 



, K ; 



k-space 



(b) 

fci(fci) 




-l> fc2(fc2) 




Fig. 2 (a) In k-space, there are two Fermi points (the K and K' points), indicated by solid circles, where the valence 
and conduction bands touch. The chiral vector specifies the lines on which energy eigenstates are located. If lines 
crossing the Fermi points (dashed arrows) exist, then the energy band gap vanishes and the system exhibits metallic 
properties. Otherwise, the system is semiconducting (solid arrows), (b) Illustration of low energy modes. The low 
energy electrons around the K-point (K'-point) generally rotate around the axis in the clockwise (counterclockwise) 
direction. The lines for armchair {Ch ~ {n, n)) nanotubes correspond to the horizontal axis in k-space; therefore, the 
low energy electrons are always classified as non-rotating states. 



3 Time-reversal symmetric gauge field In the previous section we explained why nanotubes can show 
either metallic or semiconducting behavior, one of the most remarkable properties of nanotubes. This 
is due to (1) the periodic boundary condition around the axis and (2) the energy dispersion relation of 
graphene (there are two Fermi points). Because the periodic boundary condition around the axis strongly 
affects the energy band structure, it would be interesting to consider how the solid state properties of a 
nanotube depend on the boundary condition along the axis. If the two ends of a nanotube are connected, 
the boundary condition along the axis becomes periodic. The resultant geometry, a torus, possesses an 
important structural degree of freedom called "twist" [9]. Twist can mix the motions along and around the 
axis, and is thus expected to affect the solid state properties. In Section 3.1, the kinematics of persistent 
current in a twisted carbon torus are reported, and it is shown that current exhibits the special character 
of tubule geometry. It is noted that persistent currents in a twisted torus can be well described by the 
notion of a twist-induced gauge field [4]. Except for the kinematics, the effect of a lattice deformation on 
the low-energy modes is crucial for understanding nanotube physics. In Section 3.2, the effect of lattice 
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deformation on the energy gap is examined in terms of a deformation-induced gauge field [7]. Throughout 
this section, the strong interplay between (1) and (2) has been emphasized, and it has been remarked that 
both twist- and deformation-induced gauge fields maintain the time-reversal symmetry of the system. 

3.1 Twist-induced gauge field In this section, persistent currents in toroidal carbon nanotubes are ex- 
amined. Because the effect of twist on persistent currents is of interest, the difference between persistent 
currents in an untwisted torus and that in a twisted torus is important. Moreover, it is not necessary to 
consider all of the electrons in the system. An examination of low-energy electrons is sufficient because 
the persistent currents relate to the energy spectrum flow at the Fermi level. Here, the lattice structure of 
an untwisted carbon nanotorus is specified initially, then a twisted nanotorus is considered. 

The kinematics for an untwisted torus are considered first. The wave vector of the low energy modes 
is given by k^^ = i/zjfci + H2k2, where -h/x? (— /x?) indicates clockwise (counterclockwise) rotation 
around the axis, as shown in Fig. 2(b). Because of the periodic boundary condition along the axis, we have 
'^low ' ^ = 27r/Z2. Next, the kinematics for a twisted carbon torus are considered. A specific type of torus, 
obtained by twisting an untwisted torus (A-type [4]), is used for simplicity. To express a twist degree of 
freedom, the translational vector is replaced as T — » = T + ^C/j, where t and d are integers. Here, 
t denotes the twist and d is the greatest common divisor of the two integers n and m. The wave vector 
of the low energy modes in the twisted torus is denoted by k^^ and should satisfy the foUowing periodic 
boundary conditions: 

fc±„-T^ = 27rM2, fcit • Cft = ±27:^;. (4) 
To understand the effect of twist, we consider the iimer product between and T as 

fcit • T = 277/22 - ■Ch = 2i, [n2 T . (5) 

The twist shifts the wave vector in k-space and depends on the electron motion around the axis specified by 
zb/i". It is important to note that the effect of twist can be provided by a gauge field for an untwisted torus. 
Therefore, the kinematics of a twisted torus are equivalent to that of an untwisted torus in the presence of a 
twist-induced gauge field, A^^'^*, satisfying A*^'"* • = and A*^^"* • T = 27r^/u?. It is also important 
to mention that the twist-induced gauge field does not break the time-reversal symmetry, contrary to the 
electromagnetic gauge field A*^™. This is because shifts the wave vector in one direction and depends 
not on electron motion but on electron charge e. The electromagnetic gauge field leads to the gauge 
coupUng given by 

where $1™ = -^^^ ■ T is the magnetic flux penetrating the torus and $o is the fundamental unit of magnetic 
flux. 

Just as persistent current is induced by magnetic flux $|'" or A"^, the twist-induced gauge field A}"^^^^ 
also affects persistent current. Sawtooth curves for the K-point and K' -point persistent currents are shifted 
in different directions by A*^"** because of the ± sign in Eq. (5). Suppose the persistent current in an 
untwisted torus is given by the usual sawtooth curve as a function of $1™. In this case, persistent currents 
of the K-point and K' -point are equal, and the amplitude of total current is twice that of the K-point only 
(Fig. 3(a)). The twist shifts the curves in the positive and negative directions depending on electron motion 
around the axis as shown in Fig. 3(b). The total persistent current is the sum of these two currents, giving 
a generally compUcated resulting current. Fig. 3(c) depicts a special case where the period of the current 
becomes $o/2. 
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Fig. 3 (a) Sawtooth curve of a persistent current (/pc) in an untwisted (m-class [4]) torus as a function of magnetic 
flux. The period of the curve is given by the flux quantum $o- (b) The twist shifts curves in the positive and negative 
directions depending on the clockwise and counterclockwise electron motion around the axis, (c) The total persistent 
current is given by the sum of the two currents. 

3.2 Deformation-induced gauge field We consider the effect of a lattice deformation on the low en- 
ergy modes. The Hamiltonian of an undeformed tube can be decomposed into two parts at low energy 
as Tio Hk ffi 'Hk', where Hk = vpcr ■ p (Hk' = vpa' ■ p) defines the effective dynamics for the 
K-point (K' -point) and the coiTesponding Schrodinger equation is given by the Weyl equation [10]. Here, 
the Fermi velocity is denoted by vf{^ 3V"7racc/2), Ci are the Pauli matrices and a ■ p = aipi + (T2P2 
(<t' ■ p = — (Tipi + (T2P2)- It is noted that p is the momentum measured from the K-point or K'-point. A 
lattice deformation can be included by adding the deformed Hamiltonian Tidofoim to Ho as 

■Wdoform = ^ ^SVa{ri)al_^,^a^ + h.c, (7) 

a=l,2,3 ieA 

where SVa{ri) is a modulation of the hopping amplitude from to r^+a, induced by the deformation. 
The leading order of the total Hamiltonian, Ho + TYdoform, can be expressed at low energy by Hk = 
vpa -{p ~ A) and Hk' ~ vpc' ■ {p + A), where <j ■ A — aiAi + (T2/I2, A is the vector field defined on the 
surface and vpA{r) is given by a linear function of SVa{r). It is noted that the deformation can be included 
as a gauge field because of the gauge coupling p ^ p ± A. Here, A is called the deformation-induced 
gauge field, and it does not break the time-reversal symmetry because the sign in front of A is different for 
the K and K' points. It is also important to note that this sign characteristic is equivalent to A*™"'*, but in 
this case A can be position dependent. 

The effect of A on the energy gap is now considered. For this purpose, Hk is considered in the presence 
of the electromagnetic gauge field [11, 12]. The Hamiltonian of the system is given by Hk ~ vpcr ■ 
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{p — A — eA^™), which indicates that the low energy electrons around the K-point are affected by the 
gauge field of ^ + eA^"^. The magnetic flux parallel to the tube axis is given by the sum of a (true) 
magnetic flux = § dxi A"" and a magnetic flux produced by the deformation-induced gauge field 
$deform(^2) = / dxiAi{xi,X2) as ^i{x2) = $f ^°™(a;2) + ^f™, shifting the wave vector through the 
AB effect. This shift corresponds to the fact that the lines in Fig. 2(a) also shift, resulting in a change 
of the energy gap [13, 14]. The important point here is that the change depends on the axis position X2, 
because the flux is a function of X2. Therefore, the energy gap modulates along the axis. The energy gap 
modulation can be estimated by 

It is important to note that a magnetic field appears when the energy gap is locally modulated. This 
magnetic field, perpendicular to the tube surface, can be defined by = CijdiAj = V^^^b. In fact, we 
have ^i{x2 + ^2) — ^1(2:2) = 5x2 § dxid2Ai = —SX2B3. A deformed (by Ceo) nanotube and the 
corresponding local modulation of energy gap are illustrated in Fig. 4. 

(a) A nanotube is deformed by Ceo (b) Deformation-induced gauge field A 




(c) Deformation-induced magnetic field (d) Energy gap 




Fig. 4 A deformed carbon nanotube (a) can be thought of as an undeformed nanotube in the presence of a deformation- 
induced gauge field (b). The relationship between (c) a deformation-induced magnetic field B3 (d) the corresponding 
local energy gap. It is noted that this is not the usual electromagnetic field (3>i™ = 0). 

It has been assumed thus far that the deformation does not mix the K and K' points. A short-distance 
lattice deformation, or topological defects such as a pentagon or heptagon, can mix them [15] such that 
the deformation-induced gauge field should be generahzed. A non-Abehan extension of the A field can 
successfully include these kinds of deformations [7]. 

4 Discussion The following sections discuss several problems not yet solved but which seem to be 
interesting subjects. 

4. 1 Persistent Currents in the Torus An interesting theoretical issue is the need to understand persis- 
tent currents in toroidal carbon nanotubes in more detail. There are several theoretical models for toroidal 
shapes composed of carbon atoms, classified into three groups: (tl) a small torus consisting of about hun- 
dred carbon atoms [16], (t2) a torus with a large aspect ratio in the presence of pentagon-heptagon pairs 
in the tube surface, and (t3) a torus with a large aspect ratio in the absence of pentagon-heptagon pairs. A 
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simple (t3) toras is considered here. However, even for (t3), a theoretical study of persistent currents in the 
presence of a lattice deformation is not complete. Another interesting issue is clarification of the effects of 
topological defects on persistent currents (t2). 

4.2 Local Energy Gap at the Edge This section discusses the applicability of the local energy gap 
formula. Because it is meant for non-mixing cases, it cannot be used in its present form to estimate the 
local energy gap near topological defects. For example, there may be a closed cap structure containing 
pentagons at the edge of a tube for which the local energy gap should be considered in the framework of 
the non-Abelian deformation-induced gauge field [7]. In the absence of topological defects, bonds along 
the axis direction cut at the edge, which generally gives a local deformation-induced gauge field around 
the axis Ai (see Fig. 5). Hence, it is naturally expected to cause local modulation of the energy gap at the 
edge. It should be noted that the framework for the local energy gap presented here does not give the wave 
function of energy eigenstates in a deformed nanotube. This should be clarified in future work. 



Fig. 5 Illustration of an edge part of zigzag nanotube, Ch = {n, 0). The hopping integral along the tube axis should 
be modulated at the boundary, giving rise to the deformation-induced gauge field around the axis. It is noted that a 
localized edge state exists around the zigzag edge whose energy eigenvalue is flat at the Fermi energy [17]. 
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